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INTRODUCTION 


Breast  tomosynthesis  is  a  tomographic  imaging  technique,  and  it  has  the  potential  advan¬ 
tage  to  overcome  a  major  limitation  of  conventional  mammography  through  recovering,  to  a 
large  degree,  the  loss  of  3D  information  about  the  breast  in  conventional,  2D  mammography. 
In  the  last  several  years,  there  has  been  renewed  interest  in  developing  breast  tomosynthesis 
for  detection  of  breast  cancer  [1,2].  Although  considerable  progress  has  been  made,  im¬ 
provements  to  several  areas  of  breast  tomosynthesis  technology  are  still  needed  before  it 
becomes  suitable  for  routine  clinical  use.  In  essence,  breast  tomosynthesis  can  be  considered 
as  a  dedicated  computed  tomography  with  limited  view  for  breast  imaging,  and  it  thus  re¬ 
quires  the  development  of  special  reconstruction  algorithms  for  recovering  3D  breast  images 
from  tomosynthesis  data.  In  addition,  various  physical  factors  in  breast  tomosynthesis  can 
strongly  affect  the  resulting  image  quality,  and  the  issue  of  patient  radiation  dose  in  breast 
tomosynthesis  is  of  a  concern.  The  overall  objective  of  this  project  is  to  investigate  and 
develop  reconstruction  algorithms  for  obtaining  breast  images  of  practical  use,  to  investigate 
and  evaluate  systematically  the  effects  of  various  physical  factors  on  image  quality  in  breast 
tomosynthesis,  and  to  use  and  evaluate  (empirical)  techniques  for  effectively  compensating 
for  the  effects  on  breast  tomosynthesis  images  and  for  possibly  reducing  imaging  radiation 
dose  in  breast  tomosynthesis.  It  is  fully  expected  that  the  research  will  contribute  to  the 
effort  in  the  field  to  develop  and  improve  breast  tomosynthesis  for  its  clinical  use.  This 
report  summarizes  the  progress  of  this  Predoctoral  Traineeship  Award  project  made  by  the 
recipient  during  the  past  one  year. 


BODY 


1  Training  Accomplishments 

At  the  time  of  this  report,  the  recipient,  Dan  Xia,  of  the  Predoctoral  Traineeship  Award  has 
taken  22  out  of  22  required  courses  towards  his  Ph.D.  degree  in  medical  physics.  The  courses 
include  physics  of  medical  imaging,  physics  of  radiation  therapy,  mathematics  for  medical 
physicists,  image  processing,  statistics,  anatomy  of  the  body,  radiation  biology  and  teaching 
assistant  training. 

2  Research  Accomplishments 

2.1  Investigation  of  reconstruction  algorithms  for  breast  tomosynthesis 

In  tomosynthesis,  data  are  acquired  only  at  a  small  number  of  projection  views  over  a  limited 
angular  range.  Therefore,  tomosynthesis  data  are  highly  sparse  as  compared  to  data  acquired 
in  conventional  computed  tomography  (CT).  Consequently,  existing  analytical  algorithms  for 
accurate  reconstruction  of  CT  images  are  generally  not  suitable  for  yielding  useful  images 
from  tomosynthesis  data.  For  example,  the  effect  of  filtering  may  not  completely  be  canceled 
out  due  to  the  limited  number  of  views  and  limited  angular  range,  resulting  in  prominent 
artifacts  in  reconstructed  images.  In  contrast,  iterative  algorithms  can  generally  produce 
images  with  less  artifacts  than  can  the  analytic  algorithms.  In  the  project,  we  are  inves¬ 
tigating  and  developing  a  total  variation  (TV)  based  iterative  algorithm  for  reconstructing 
accurate  images  from  incomplete  projection  data  [3].  We  have  investigated  the  following 
issues:  (1)  Performance  of  the  TV-based  algorithm  when  applied  to  tomosynthesis  data  gen¬ 
erated  from  discrete  images.  In  this  case,  one  has  a  discrete  matrix  system,  which  allows 
ideal  reconstructions  of  the  underlying  discrete  images.  The  result  in  this  study  provides  the 
upper  bound  on  the  performance  of  the  TV-based  and  any  other  reconstruction  algorithm. 
(2)  Convergence  conditions  of  the  TV-based  algorithm  under  different  data  conditions  and 
different  constraint  parameters.  This  investigation  is  critical  for  any  iterative  algorithm  to 
be  practically  useful.  We  have  successfully  addressed  issue  (1)  and  are  working  actively  on 
issue  (2). 

2.2  Investigation  of  scanning  configurations  in  breast  tomosynthesis 

Currently,  breast  tomosynthesis  acquires  data  at  about  12  to  20  projection  views  over  a 
limited  angular  range  around  20  to  50  degrees.  Typically,  a  circular  source  trajectory  is 
adopted  for  collection  of  cone-beam  projections.  It  remains,  however,  largely  unexplored 
as  to  what  the  optimal  scanning  configurations  are.  In  the  last  several  months,  we  have 
conducted  initial  investigation  on  image  reconstruction  from  data  acquired  (1)  at  a  small 
number  of  views  and  (2)  over  a  limited  angular  range.  In  the  studies,  data  were  collected 
from  different  phantoms,  including  a  breast  phantom  developed  in  Ref.  [4],  at  different 
numbers  of  views. 

We  first  generated  cone-beam  data  at  15,  20,  40,  and  60  projection  views  over  27r.  The 
reason  for  the  angular  range  is  27r  is  that  the  study  results  would  not  be  affected  by  the  issue 
of  the  limited  angular  range.  From  these  data  sets,  we  have  used  the  TV-based  algorithm 
to  reconstruct  images.  For  comparison,  we  have  also  conducted  preliminary  reconstructions 
by  using  the  expectation-maximization  (EM)  algorithm  and  the  algebraic  reconstruction 


technique  (ART).  Although  we  are  currently  in  the  process  of  analyzing  these  reconstruc¬ 
tion  results,  the  initial  evaluation  appears  to  indicate  that  the  TV-based  algorithm  is  more 
accurate  than  the  EM  and  ART  algorithms  in  reconstructing  images  from  few-view  data. 

We  subsequently  studied  image  reconstruction  from  data  acquired  over  a  limited  angular 
range.  Specifically,  we  have  selected  angular  ranges  of  7t/4,  7t/2,  and  n  and  collected  data 
at  15,  20,  40,  and  60  views  over  each  of  these  angular  ranges.  From  the  collected  data, 
we  reconstructed  images  by  use  of  the  TV-based  algorithm.  For  comparison,  the  EM  and 
ART  algorithms  were  also  used  for  reconstructing  images  from  these  data  sets.  Our  initial 
analysis  of  these  results  suggests  [5]  that  (1)  when  the  angular  range  decreases,  image  quality 
obtained  with  these  algorithms  decreases,  and  (2)  the  TV-based  algorithm  generally  yields 
images  with  less  artifacts  than  do  the  EM  and  ART  algorithms.  We  are  still  in  the  process 
of  conducting  additional  studies  on  this  problem. 

2.3  Investigation  of  the  physical  factors  in  tomosynthesis  imaging 

Various  physical  factors  can  significantly  affect  image  quality  in  breast  tomosynthesis.  These 
factors  include  data  noise,  non-uniform  image  spatial  resolution,  scatter,  and  detector  re¬ 
sponses.  We  have  begun  to  perform  research  on  investigating  and  correcting  for  the  effect  of 
physical  factors  on  image  quality.  In  our  preliminary  studies,  we  have  evaluated  the  effects 
of  two  physical  factors,  (1)  cone-beam  data  noise  and  (2)  non-uniform  image  resolution. 

In  an  attempt  to  separate  the  issues  of  few-view  and  limited  angular  range  from  data  noise, 
we  have  first  performed  a  thorough  noise  study  for  image  reconstruction  from  projection 
data  collected  at  a  large  number  of  views  over  2ti.  The  result  of  this  study  will  provide  a 
theoretical  guidance  to  the  investigation  of  the  noise  properties  in  tomosynthesis  in  which 
data  are  acquired  at  a  small  number  of  views  over  a  limited  angular  range.  In  this  study, 
we  specifically  focus  on  investigating  how  data  noise  are  propagated  into  the  reconstructed 
images.  Both  analytic  and  numerical  analyses  were  carried  out,  and  the  results  indicate  that 
variances  of  reconstructed  images  are  spatially  varying  and  that  the  levels  of  variances  in 
different  regions  are  not  affecting  significantly  each  other  [6].  Because  the  study  was  based 
upon  analytic  result,  it  was  computationally  possible  to  accomplish  this  initial  study.  On 
the  other  hand,  noise  studies  involve  the  iterative  algorithms  such  as  the  TV-based,  EM, 
and  ART  algorithms  are  much  more  demanding  computationally.  We  have  just  begun  to 
perform  the  image  noise  studies  by  using  the  iterative  algorithms. 

In  current  breast  tomosynthesis,  the  spatial  resolution  within  a  transverse  plane  is  much 
finer  than  that  along  the  longitudinal  direction,  and  image  representation  with  non-isotropic 
spatial  resolution  is  used  in  iterative  algorithms  for  reducing  computational  time.  Such  an 
image  representation  can  lead  to  significant  artifacts  in  iterative  reconstruction.  In  the  last 
couple  of  months,  we  have  been  investigating  the  effect  of  non-isotropic  image  representation 
on  iterative  reconstruction  accuracy  of  breast  tomosynthesis  images.  We  have  reconstructed 
images  by  use  of  TV-based,  EM,  and  ART  algorithms  for  image  representations  with  different 
ratios  of  the  in-plane  and  longitudinal  resolution.  Our  results  demonstrate  that  non-isotropic 
image  representation  can  lead  to  significant  artifacts  in  reconstructed  images.  The  appear¬ 
ance  and  severity  of  the  artifacts  depend  not  only  upon  the  ratio  between  the  in-plane  and 
longitudinal  resolution  but  also  upon  the  iterative  algorithms.  The  TV-based  algorithm 
seems  to  be  less  susceptible  to  the  effect  than  the  EM  and  ART  algorithms.  Through  the 
selection  of  algorithm  parameters,  the  artifacts  can  be  reduced  [7]. 


KEY  RESEARCH  ACCOMPLISHMENTS 


We  have  implemented  the  modification  to  the  TV-based  algorithm  tailored  to  breast 
tomosynthesis. 

We  have  also  implemented  and  tested  the  EM  and  ART  algorithms. 

We  have  conducted  an  investigation  of  the  upper  bound  on  the  performance  of  TV-based 
algorithm  for  image  reconstruction  from  tomosynthesis  data. 

We  have  conducted  a  preliminary  investigation  of  the  convergence  property  of  the  TV- 
based  algorithm  under  different  data  conditions  and  different  constraint  parameters  in 
reconstruction  of  tomosynthesis  images. 

We  have  carried  out  preliminary  investigations  of  the  scanning  configurations  in  breast 
tomosynthesis  imaging  by  using  different  numbers  of  views  over  different  angular  ranges. 

We  have  performed  a  preliminary  study  on  the  effect  of  some  physical  factors  on  breast 
tomosynthesis.  The  physical  factors  were  considered  in  the  study  include  data  noise 
and  non-isotropic  spatial  resolution. 
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CONCLUSIONS 


The  recipient  of  the  Predoctoral  Traineeship  Award  has  finished  the  required  courses 
towards  his  Ph.D.  degree.  These  trainings  have  proven  useful  for  the  recipient  to  achieve  the 
proposed  research  goals. 

During  the  first  year,  we  have  implemented  and  investigated  the  TV-based  algorithm  for 
image  reconstruction  in  breast  tomosynthesis.  We  first  evaluated  the  performance  of  the 
algorithm  under  ideal  conditions  and  its  convergence  properties.  We  have  subsequently  con¬ 
ducted  numerical  studies  to  investigate  the  image  reconstruction  by  use  of  the  TV-based 
algorithm  and  the  existing  EM  and  ART  algorithms  for  different  imaging  configurations  in 
breast  tomosynthesis.  Our  results  suggest  that,  in  general,  the  TV-based  algorithm  pro¬ 
vide  tomosynthesis  images  with  quality  higher  than  EM  and  ART  algorithms  in  terms  of 
the  metrics  such  as  root-mean-square  error  and  image  resolution.  Furthermore,  we  have 
investigated  the  effect  of  some  physical  factors,  such  as  data  noise  and  non-isotropic  spatial 
resolution,  on  tomosynthesis  images.  Overall,  we  have  achieved  the  goals  for  the  first  year 
and  laid  down  the  foundation  for  the  research  in  the  next  two  years.  Our  aims  in  the  next 
two  years  include,  based  upon  what  we  have  learned  from  the  results  described  above,  fur¬ 
ther  investigating  and  streamlining  the  TV-based  algorithm  in  terms  of  convergence,  speed, 
and  rubostness.  Careful,  comparative  evaluation  studies  of  reconstruction  algorithms  will 
be  conducted  under  more  realistic  physical  conditions.  In  particular,  we  will  expand  our 
data-generation  programs  to  include  non-linear  partial  volume,  beam-hardening,  and  scatter 
effects  in  our  data.  We  will  implement  existing  algorithms  to  correct  for  the  scatter  effect 
before  image  reconstruction.  Additionally,  we  will  generalize  reconstruction  algorithms,  in¬ 
cluding  the  TV-based,  EM,  and  ART  algorithms  to  incorporate  these  factors  so  that  they 
can  be  compensated  for  during  reconstruction  iterations.  We  will  look  into  the  results  from 
current  studies  in  other  groups  in  our  department  on  the  effect  of  detector  response,  which 
is  expect  to  be  spatially  varying,  on  data  acquisition.  If  it  is  significant,  we  will  include  the 
detector  response  into  the  data-generation  program  and  reconstruction  algorithms.  We  will 
also  perform  tomosynthesis  experiments  to  collect  real  data  of  physical  breast  phantoms  for 
additional  evaluation  of  the  scanning  configurations  and  reconstruction  algorithms.  Finally, 
we  will  select  real-patient  data  for  testing  and  assessing  the  reconstruction  algorithm. 
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Clinical  Relevance/ Application: 

Breast  tomosynthesis  has  received  renewed  interest  because  it  can  provide  3D 
information  about  the  breast.  This  work  concerns  iterative  reconstruction  of  accurate 
breast  tomosynthesis  images. 

Purpose: 

In  current  breast  tomosynthesis,  image  representation  with  non-isotropic  spatial 
resolution  is  used  for  reducing  computational  time.  This  can,  however,  lead  to  artifacts  in 
iterative  reconstruction  of  breast  tomosynthesis  images.  In  the  work,  we  investigate  the 
effect  of  non-isotropic  image  representation  on  the  reconstruction  accuracy.  Based  upon 
the  investigation,  we  devise  schemes  for  reducing  artifacts  in  iterative  reconstruction. 

Materials  and  Method: 

In  the  work,  we  focus  on  investigating  the  effect  of  non-isotropic  image  representation  on 
reconstruction  accuracy  of  iterative  algorithms.  The  iterative  algorithms  under  study 
include  the  total- variation  (TV)  based,  expectation  maximization  (EM),  and  algebraic 
reconstruction  technique  (ART)  algorithms.  Tomosynthesis  data  are  generated  at  12  and 
20  views  over  50  degrees  from  phantoms,  including  a  breast  phantom.  We  have 
reconstructed  images  by  using  image  representations  with  different  degrees  of  non¬ 
isotropic  spatial  resolution.  Specifically,  in  each  image  representation,  the  ratio  between 
the  in-plane  and  longitudinal  resolution  for  an  image  voxel  is  selected  to  be  a  value  less 
than  1 . 

Results: 

We  have  reconstructed  images  by  use  of  TV-based,  EM,  and  ART  algorithms  for  image 
representations  with  different  ratios  of  in-plane  and  longitudinal  resolution.  Our  results 
demonstrate  that  non-isotropic  image  representation  can  lead  to  significant  artifacts  in 
reconstructed  images.  The  appearance  and  severity  of  the  artifacts  depend  not  only  upon 
the  ratio  between  the  in-plane  and  longitudinal  resolution  but  also  upon  the  iterative 
algorithms.  The  TV-based  algorithm  seems  to  be  less  susceptible  to  the  effect  than  the 
EM  and  ART  algorithms.  Through  the  selection  of  algorithm  parameters,  the  artifacts  can 
be  reduced. 

Conclusion: 

The  non-isotropic  image  representation  can  significantly  affect  reconstruction  accuracy 
obtained  with  iterative  algorithms  in  breast  tomosynthesis. 
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Purpose: 

Investigation  and  application  of  scanning  and  reconstruction  approaches  to  tomosynthesis 
breast  imaging. 

Materials  and  Method: 

Tomographic  imaging  such  as  tomosynthesis  has  the  advantages  to  overcome  the  major 
limitation  of  conventional  mammography.  The  loss  of  information  in  2D  projection 
imaging  can  be  recovered  in  tomographic  imaging.  In  tomosynthesis,  the  acquired  data 
are  highly  sparse  as  compared  to  CT  data.  Therefore,  image  reconstruction  in 
tomosynthesis  is  challenging.  Moreover,  it  remains  largely  unexplored  as  to  what  the 
optimal  scanning  configurations  are.  Recently,  we  have  developed  a  TV-based  iterative 
algorithm  for  accurate  images  from  sparse  or  incomplete  projection  data.  In  this  work,  we 
will  investigate  image  reconstruction  by  use  of  the  TV-based  algorithm.  Moreover,  using 
this  algorithm  and  other  algorithms,  we  investigate  the  optimal  scanning  geometries  in 
tomosynthesis  breast  imaging.  Specifically,  we  have  investigated  scanning  configurations 
in  terms  of  number  of  projection  views  and  the  range  of  scanning  angle. 

Results: 

We  have  conducted  numerical  studies  to  investigate  image  reconstruction  and  scanning 
parameters  by  use  of  the  TV-based  algorithm  and  other  algorithms.  From  numerical  and 
physical  phantom,  we  have  acquired  projection  data  from  physical  phantoms  for  different 
scanning  configurations.  From  these  data  sets,  we  reconstruct  images  by  using  the  TV- 
based  algorithm.  For  comparison,  other  algorithms,  including  the  FBP  algorithm.  Using 
various  evaluation  indices,  we  characterize  the  reconstructed  image  quality  and,  from 
which,  we  evaluate  the  performance  of  the  reconstruction  algorithms  and  scanning 
configurations. 

Conclusion: 

In  this  work,  we  have  investigated  different  scanning  configurations  and  reconstruction 
algorithms,  including  the  TV-based  algorithm,  in  tomosynthesis  breast  imaging.  Both 
computer-simulation  data  and  real  physical  phantom  data  were  used.  Our  results  suggest 
that,  in  general,  the  TV-based  algorithm  can  provide  higher  image  quality  for 
tomosynthesis  imaging  than  other  algorithms  and  that  increasing  the  scanning  angular 
range  can  considerably  improve  image  quality. 
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Abstract 

Recently,  there  has  been  much  progress  in  algorithm  development  for  image  reconstruction  in  cone-beam 
computed  tomography  (CT).  Current  algorithms,  including  the  chord-based  algorithms,  now  accept  minimal  data 
sets  for  obtaining  images  on  volume  regions-of-interest  (ROIs)  thereby  potentially  allowing  for  reduction  of  x- 
ray  dose  in  diagnostic  CT.  As  these  developments  are  relatively  new,  little  effort  has  been  directed  at  investigating 
the  response  of  the  resulting  algorithm  implementations  to  physical  factors  such  as  data  noise.  In  this  work,  we 
perform  an  investigation  on  the  noise  properties  of  ROI  images  reconstructed  by  using  chord-based  algorithms 
for  different  scanning  configurations.  We  find  that,  for  the  cases  under  study,  the  chord-based  algorithms  yield 
images  with  comparable  quality.  Additionally,  it  is  observed  that,  in  many  situations,  large  data  sets  contain 
extraneous  data  that  may  not  reduce  the  ROI-image  variances. 
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1  Introduction 

In  recent  years,  exact  algorithms  have  been  developed  for  reconstructing  images  [1]  and  for  reconstructing  images 
on  “7r-lines”  [2,  3,  4]  from  helical  cone-beam  data.  Since  2005,  works  have  being  published  on  algorithm  devel¬ 
opment  for  reconstructing  images  on  chords  for  general  trajectories  [5,  6,  7,  8],  Some  of  these  algorithms  can 
reconstruct  images  within  three  dimensional  (3D)  regions  of  interest  (ROIs)  from  cone-beam  data  containing  both 
longitudinal  and  transverse  truncations.  The  introduction  of  the  A/ -line  concept  and  reconstruction  [5,  9]  provides 
additional  flexibility  for  covering  volume  ROIs. 

As  these  algorithm  developments  are  relatively  recent,  little  effort  has  been  directed  at  investigating  their  noise 
properties.  With  the  algorithm  development  for  ROI-image  reconstruction,  it  has  been  tacitly  assumed  that  the 
reduction  in  necessary  scanning  angle  and  in  projection  data  may  lead  to  ROI  images  from  less  radiation  exposure. 
This  conclusion  may,  however,  depend  on  the  noise  properties  of  reconstruction  algorithms.  If  ROI  reconstruction 
from  the  minimal  (or  reduced)  data  set  leads  to  noisier  ROI  images  than  reconstruction  of  the  same  ROI  from  a 
larger  data  set,  it  may  be  necessary  to  increase  the  x-ray  source  intensity  for  the  ROI-data  set  to  attain  the  same 
image  quality  as  those  reconstructed  from  larger  data  sets.  Such  an  increase  can  offset  the  fact  that  reduced  or 
minimum  projection  data  are  needed  for  ROI  reconstruction. 

The  focus  of  this  work  is  to  investigate  the  noise  properties  of  image  reconstruction  from  minimal  data  set 
and  large  data  sets  by  use  of  chord-based  algorithms.  We  demonstrate  that  the  minimal  data  set  can  indeed 
lead  to  actual  reduction  of  radiation  exposure  for  attaining  comparable  image  quality,  defined  in  terms  of  image 
variance,  as  that  obtained  with  a  larger  data  set.  In  Sec.  2,  we  briefly  summarize  the  chord-based  reconstruction 
algorithms:  backprojection  filtration  (BPF)  [2,  6],  minimum  data  filtered  backprojection  (MDFBP)  [4,  6],  and 
filtered  backprojection  (FBP)  [6,  10]  algorithms.  In  Sec.  3,  we  perform  analysis  and  empirical  studies  on  noise 


properties  of  images  reconstructed  from  parallel-beam  fan-beam,  and  cone -beam  data.  Finally,  discussion  is  given 
in  Sec.  4. 


2  Chord-based  Reconstruction  Algorithms 


We  consider  a  continuous  source  trajectory  specified  by  ro(s)  =  ( x(s ),  y(s),  z(s)),  where  x(s),  y(s),  and  z(s) 
denote  the  x-,  y-,  and  ^-components  of  ro(s)  in  the  fixed-coordinate  system,  and  s  is  a  curve  parameter  indi¬ 
cating  the  position  of  the  x-ray  source  on  the  trajectory.  The  projection  data  of  the  object  function  f(r)  can  be 
mathematically  expressed  as 
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D(r0(s),P)  =  /  dt  f(r0{s) +t(. 3),  (1) 

Jo 


where  the  unit  vector  j3  denotes  the  direction  of  a  specific  x-ray  passing  through  the  point  r.  We  also  introduce  two 
additional  coordinate  systems  {it,  v.  u>}  and  {ud,  vd}  to  describe  the  geometry  in  a  general  scan.  They  are  fixed 
on  the  rotating  source  point  and  the  cone-beam  projection  of  the  source  point,  respectively,  which  are  referred 
to  as  the  rotation-coordinate  and  detector-coordinate  systems.  Let  eu(s),  e„(s),  and  ew(s)  denote  the  orthogonal 
unit  vectors  of  the  rotation-coordinate  system.  The  rotation-coordinate  system  can  be  chosen  such  that  eu(s) 
and  ew(s)  are  within  the  x-y  plane  and  ev(s)  is  parallel  to  the  .---axis.  One  can  also  choose  the  “well  oriented" 
coordinate  system  as  the  rotation-coordinate  system  [5],  in  which  unit  vector  eu(s)  is  parallel  to  and  unit  vectors 
ev(s)  and  ew(s)  are  perpendicular  to  the  direction  of  dr°^ .  We  assume  that  a  detector  plane  is  placed  at  a 
distance  S  from  the  source  point  and  orients  along  ew(s).  The  detector-coordinate  system  {ud,  Vd}  is  the  cone- 
beam  projection  of  the  2D-coordinate  system  {it,  v }  onto  the  detector  plane,  and  the  Ud~  and  tv/-axis  are  along 
e„(s)  and  ev(s),  respectively.  In  this  situation,  we  also  use  P(ud,  Vd,  s)  to  denote  the  cone-beam  projection,  thus 
D(r0(s),j3)  =  P(ud,Vd,s),  when 
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[udeu(s)  +  vdev(s)  -  Sew(s))  and  A(ud,  vd)  =  \f  u2d  +  vj  +  S2. 
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In  a  2D  case,  it  can  be  observed  that  vd  =  0.  For  notational  convenience,  we  use  A(ud)  and  P(ud,  s)  to  denote 
A(ud,  0)  and  P(ud ,  0,  s),  respectively. 

A  chord  is  a  line  segment  connecting  two  points  ro(sa)  and  ro(s(,)  on  the  trajectory.  Any  point  r  on  the  chord 
can  be  expressed  as 


ro{sa)  +  r0{sb) 

- ^ - 1-  xc  e, 
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where  ec  =  |  (s*1 )  |  ^enotes  the  direction  of  the  chord,  and  l  =  2  |ro(sb)  —  ro(sa)|  is  one  half  of  the  chord 

length.  For  a  helical  trajectory,  the  curve  parameter  s  is  linearly  related  to  the  rotation  angle  A,  and  in  the  current 
work,  we  select  s  =  A.  When  sa  and  st,  are  within  one  turn,  the  chord  becomes  the  conventional  7r-line  segment 
[2,  11,  12],  The  intersection  between  a  chord  and  the  object  is  referred  as  a  support  segment.  Let  xc\  and  xC2 
represent  the  end  points  of  a  support  segment.  Because  the  trajectory  under  consideration  never  intersects  the 
object,  we  have  [xci,xc2\  C  [—1,  l ].  Therefore,  one  can  use  (xc,  sa,  Sb )  and  fc{ xc,  sa,  Sb )  to  denote  a  point  and 
the  corresponding  image  on  the  chord.  We  have  previously  developed  three  algorithms,  which  are  referred  to  as 
the  BPF  [2,  6,  10],  MDFBP  [4,  6],  and  FBP  [6,  10]  algorithms,  respectively,  for  exact  image  reconstruction  on  a 
chord  of  a  general  trajectory. 


2.1  The  BPF  algorithm 

The  BPF  algorithm  [2,  6]  reconstructs  the  image  on  a  chord  specified  by  sa  and  as 
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where  xc  €  [xa,  is],  and  parameters  xa  and  xb  are  two  points  on  the  chord  satisfying  [xci,  xC2]  C  [xa,  xR]  C 
[— l,  l ].  The  function  6(xc)  is  defined  as  b(xc)  =  -\/ (xb  —  xc)(xc  —  xa),  and  P(udo ,  vdo,  sa)  denotes  the  projec¬ 
tion  along  the  chord  specified  by  sa  and  The  filtered  image  /(xc,  sa,  sb)  is  given  by 
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where  the  backprojection  image  on  the  chord  can  be  written  as 
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and  the  rect  function  IIc(x' )=1  if  xj,  €  [xa,  Xb]  and  0  otherwise.  It  can  be  observed  in  Eq.  (4)  that  the  chord 
image  can  be  obtained  exactly  from  knowledge  of  the  backprojection  image  g(x'c,  sa,  sb)  for  x'c  £  [xa,xb], 
which  we  refer  to  as  the  reconstruction  segment  because  it  determines  the  actual  reconstruction  interval  on  the 
chord.  In  particular,  because  the  reconstruction  segment  [xa,Xb\  can  be  chosen  as  small  as  the  support  segment 
[xci,xC2],  the  chord  image  can  be  reconstructed  from  knowledge  of  g(x'c,  sa,  sb)  only  on  the  support  segment. 
This  interesting  property  of  the  Hilbert  transform  forms  the  basis  for  exact  image  reconstruction  on  a  chord  from 
projections  containing  longitudinal  or  transverse  truncations  [13]. 


2.2  The  MDFBP  algorithm 

The  BPF  algorithm  reconstructs  the  chord  image  by  performing  a  ID  filtration  (i.e.,  the  integration  over  x'c  in  Eq. 
(4))  of  the  backprojection  image  (i.e.,  the  integration  over  s  in  Eq.  (6)).  On  the  other  hand,  the  MDFBP  algorithm 
[4,  6]  reconstructs  the  chord  image  by  performing  a  ID  data  filtration  (i.e.,  the  integration  over  u'c)  prior  to  their 
backprojection  (i.e.,  the  integration  over  s)  onto  the  chord: 
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where  the  modified  data  function  is  given  by 
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w\  =  —  [r'o(so)  —  ro(s)]  •  e^,,  and  w2  =  —  [rb(s;,)  —  ro(s)]  •  ew-  For  a  source  position  s,  the  variables  uc  and  u'c 
denote  the  cone-beam  projections  of  xc  and  x'c  onto  the  detector  and  can  be  obtained,  respectively,  by  replacing  x 
with  xc  and  x'c  in 
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The  rect  function  Ac(x'c)  in  Eq.  (8)  indicates  that  the  MDFBP  algorithm  can  reconstruct  a  chord  image  from 
knowledge  of  data  only  on  the  cone-beam  projection  of  the  reconstruction  segment  \xa,xb\,  which  can  be  as 
small  as  the  support  segment.  Therefore,  similar  to  the  BPF  algorithm,  the  MDFBP  algorithm  can  also  reconstruct 
a  chord  image  from  data  containing  truncations  [4,  6], 


2.3  The  FBP  algorithm 

The  chord-based  FBP  algorithm  [6,  10]  can  be  expressed  as 

,  (9) 

q—s 

where  uc  indicates  the  cone-beam  projection  of  xc  onto  the  detector  and  is  determined  by  using  xc  to  replace 
x  in  Eq.  (8),  and  A  denotes  the  distance  from  the  source  point  rb(s)  to  a  point  on  the  detector  at  which  the 
ray  connecting  r  and  fo(s)  intersects  the  detector.  As  the  filtering  (i.e.,  the  integration  over  u'c)  is  carried  out 
over  the  projection  of  the  straight  line  containing  the  chord,  similar  to  other  existing  FBP-based  algorithms,  the 
chord-based  FBP  algorithm  cannot  exactly  reconstruct  ROI  images  from  data  containing  transverse  truncations. 

2.4  Data-sufficiency  conditions 

As  shown  in  Eq.  (9),  a  data-sufficiency  condition  for  the  FBP  algorithm  is:  (a)  data  are  available  over  the  trajectory 
segment  s  £  [s„,  *•/,],  and  (b)  for  each  s,  data  on  the  cone-beam  projection  of  the  chord  are  non- truncated.  This 
condition  is  similar  to  that  for  other  FBP-based  algorithms  [1,  9,  14,  15,  16],  From  Eqs.  (4)  and  (7),  a  data- 
sufficiency  condition  for  the  chord-based  BPF  and  MDFBP  algorithms  is:  (a)  data  are  collected  over  the  trajectory 
segment  [sa,  sj,] ,  and  (b)  at  each  s,  data  only  on  the  cone-beam  projection  of  the  reconstruction  segment  [a ‘a,  Xb] 
on  the  chord  are  available.  It  follows  that,  because  the  reconstruction  segment  [xa,  Xb]  can  be  chosen  as  small  as 
the  support  segment  [xci,xc2],  the  BPF  and  MDFBP  algorithms  require,  at  each  s,  data  only  on  the  cone -beam 
projection  of  the  support  segment  [xci,xc2\  (instead  of  the  entire  chord-line  as  the  chord-based  FBP  algorithm 
requires).  Different  selections  of  the  reconstruction  segment  [xa,  Xb]  imply  that  different  amounts  of  data  at  each 
s  can  be  used  for  reconstructing  the  chord  image.  Under  the  ideal  continuous  conditions,  different  selections  of 
[xa,Xb\  yield  identical  chord  images.  However,  when  data  contain  noise  and  other  inconsistencies,  and  when 
different  selections  of  [xa,xb\  are  used,  the  BPF  and  MDFBP  algorithms  in  their  discrete  forms  may  yield 
different  chord  images.  This  is  an  issue  that  will  be  investigated  below. 
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3  Noise  Properties  of  Chord-based  Image  Reconstruction 

The  BPF,  MDFBP,  and  FBP  algorithms  described  above  can  be  applied  to  reconstructing  chord  images  from 
parallel-,  fan-,  and  cone -beam  data  [17],  Algorithms  analogous  to  the  BPF  algorithm  that  are  capable  of  recon¬ 
structing  2D  ROI  images  from  truncation  data  have  also  previously  been  proposed  [8,  13,  18].  We  study  below  the 
noise  properties  of  chord-based  reconstruction  by  use  of  these  algorithms  in  their  discrete  forms.  As  mentioned 
above,  the  BPF  and  MDFBP  algorithms  can  reconstruct  the  image  on  the  reconstruction  segment  [xa,  xb]  as  long 
as  it  covers  the  support  segment  [xci,xc2\.  We  analyze  image-noise  properties  on  reconstruction  segments  of 
different  lengths. 


3.1  Analysis  of  image-noise  properties 


The  chord-based  algorithms  invoke  three  major  mathematical  operations:  differentiation,  backprojection,  and 
filtration.  To  a  large  extent,  the  BPF,  MDFBP,  and  FBP  algorithms  differ  in  the  orders  of  invoking  these  operations. 
Below,  we  focus  on  investigating  the  noise  properties  of  differentiation,  backprojection,  and  filtration  in  the  BPF 
algorithm.  The  approach  taken  in  the  investigation  is  readily  applicable  to  analyzing  the  noise  properties  of 
the  MDFBP  and  FBP  algorithms.  In  the  presence  of  data  noise,  the  measured  projection  D(r(J  (.sj,  (3)  should 
be  interpreted  as  a  stochastic  process.  (Throughout  the  paper,  we  use  boldface  and  normal  letters  to  denote 
a  stochastic  process  and  its  mean,  respectively.)  Because  the  backprojection  g (xc,sa,Sb)  and  the  final  image 
fc(xc,  s0,  Sb)  on  a  chord  are  computed  from  D(ro(s),  0),  they  should  also  be  considered  as  stochastic  processes. 
We  focus  in  this  section  on  investigating  the  chord-image  variance  ic(xc,  sa,  Sb)  by  the  investigation  of  noise 
propagation  through  each  step  involved  in  the  BPF  reconstruction  algorithm. 


3.1.1  Noise  properties  of  the  differentiation/backprojection  for  parallel-beam  data 

Let  P(ud,  s)  denote  the  parallel-beam  projection  on  detector  bin  ad  acquired  at  view  s.  We  assume  the  covariance 
of  the  projection  data  P(ud,  s)  to  be  uncorrelated,  i.e., 

Cov{P(ud,  s),  P (u'd,  s')}  =  cr2(ud,  s)6(ud  -  u'd)S(s  -  s'),  (10) 


where  a2(ud ,  s)  denotes  the  variance  of  the  projection  data.  The  backprojection  image  on  the  chord  is  given  by 

[13] 


g(xc,  Sa ,  Sb) 


ds  - — P (ud,  s ) 
oud 


(11) 


where  ud  —  r  ■  eu(s).  The  final  image  variances  on  a  chord  depend  upon  the  covariance  of  the  backprojection 
image,  which,  using  Eq.  (11),  can  be  written  as 


Cov{g(xc,  Sa,  Sb),  g(x'c,  Sa,  Sfc)} 


*'Cov|^P(“''s)'^PW’s')}- 


(12) 


The  evaluation  of  the  backprojection-image  covariance  involves  the  data-derivative  covariance,  which  can  be 
conveniently  written  as 


d  d 

Cov{— — P (ud,  s),  (u'd,  s')}  =  au>  <r2(ud,  s)S(ud  -  u’d)5(s  -  s')  +  Tp 

oud  oud 


1  i.Ud,  Ud,  S,  S  ), 


(13) 


where  Var{P(ud,  s)}  =  a2(ud,  s)  denotes  the  known  data  variance,  which  can  be  a  function  of  ud  and  s.  The 
second  term  T para{ud,  u'd,  s,  s')  represents  the  difference  between  the  term  on  the  left-hand  side  and  the  first  term 
on  the  right-hand  side  of  Eq.  (13).  Although  the  magnitude  of  Tpara{ud,  u'd,  s,  s')  can  be  larger  than  or  comparable 
to  that  of  the  first  term  in  the  right-hand  side  of  Eq.  (13),  numerical  results  below  show  that  its  contribution  to  the 
final  image  variance  on  a  chord  is  negligibly  small.  Therefore,  we  consider  only  the  first  term  in  the  derivation  of 
the  chord-image  variance  below.  The  parameters  a  and  oj  are  introduced  to  account  for  the  interpolation  effect  of 
the  discrete  data  derivative  and  discrete  backprojection  on  the  chord-image  variance.  Using  the  first  term  in  Eq. 
(13),  we  can  rewrite  Eq.  (12)  as 


Cov{g(xc,  Sa,  Sb),  g(x'c,  Sa,  Sh)} 


auj  /  ds  a2 [ud,  s)5{ud  -  u'd). 


(14) 


We  now  consider  two  points  xc  and  x'c  on  the  chord  and  let  ud  and  u'd  denote  their  parallel-beam  or  fan-beam 


projections  onto  the  detector.  Clearly,  for  a  source  position  s  that  satisfies  s  /  s„  or  sb,  one  can  conclude  that 


Ud  —  u'd  =  0  if  xc  —  x'c  =  0 

Ud  l^d  ^  0  if  Xc  Xc  ^  0. 


Thus,  if  xc  =  x'c,  Cov{g(;rc,  s0,  sb),  g(;rc,  sa,  sfc)}  =  (5(0)  aw fj„  dsa2(ud,  s)  and,  ifa’c  ^  Cov{g(xc,  sa,  sb),g(x'c,  sa,  sb) 

0.  Therefore,  the  covariance  of  the  backprojection  image  for  parallel-beam  projections  can  be  re-expressed  as 


Cov{g(xc,  sa,  sb),  g(x'c,  sa,  s6)}  «  c( xc)  S(xc  -  x'c), 


c(xc)  =  au  /  dsa2(ud,s). 


3.1.2  Noise  properties  of  the  differentiation/backprojection  for  fan-beam  data 


In  the  fan-beam  case,  we  use  P (ud,  s)  to  denote  the  projection  on  detector  bin  Ud  acquired  at  view  s.  Again,  we 
assume  P(ud,  s )  to  be  uncorrelated  and  satisfy  Eq.  (10).  The  backprojection  image  in  Eq.  (6)  can  be  re-expressed 


rsb  i 

g(a •c,sa,sb)  =  J^  ds  _.o(g)|2P>rf,s), 


Sr-e  (s') 

where  Ud  =  is  the  fan-beam  projection  of  xc  onto  the  detector.  The  weighted-data  derivative 


P \ud,  s)  is  given  by 


P'(wd,s)  =  A2(ud) 


2  .  I  dfo(s)  I  d  \P(ud,s) 


ds  |  dud  [  A(ud)  \ 


Using  Eq.  (18),  one  can  write  the  covariance  of  the  backprojection  image  as 


psb  rsb 

Cov{g(;rc,  sa,  sb),  g{x'c,  sa,  Sb)}  =  /  ds  ds' 

J  Sn  J  S, 7 


JSa  Jsa  '  |r-r0(s)|2  |r-r0(s')|2 
x  CovJP^Md,  s),  P^it^,  s')}, 


which  depends  upon  the  covariance  of  P'  (ud,  s).  Again,  we  can  conveniently  write  the  covariance  of  P  (ud,  s)  as 


Cov{P'(«d,  s),P'(u^,  s')}  =  auj<r2(ud,s)  A2(ud)  dr^s'1  S{ud  -  u'd)S{s  -  s') 

as 


+  Tfan(ud,  Cld,  S,  S  ), 


where  Var{P(itd,  s)}  =  a2(ud,  s)  denotes  the  known  data  variance,  which  can  be  a  function  of  ud  and  s.  The 
second  term  T fan(ud ,  u'd,  s,  s')  represents  the  difference  between  the  term  on  the  left-hand  side  and  the  first  term 
in  the  right-hand  side  of  Eq.  (21).  As  numerical  results  below  indicate,  it  turns  out  that  T fan(ud ,  u'd ,  s,  s')  will 
also  have  a  negligible  contribution  to  the  chord-image  variance.  Therefore,  we  consider  only  the  first  term  in  the 
derivation  of  the  chord-image  variance  below.  Again,  the  parameters  a  and  u>  are  introduced  to  account  for  the 
interpolation  effect  of  the  discrete  data  derivative  and  discrete  backprojection  on  the  chord-image  variance.  Using 
the  first  term  in  Eq.  (21),  we  can  rewrite  Eq.  (20)  as 


r,  t  ,  ^  t  ,  M  fSb  J  0- 2{ud,s)A2{ud)  dfo(s)  , 

Cov{g{xc,sa,sb),g(xc,sa,sb)\  «  auj  ds — F _  ,  ,, - - —  d(ud-ud).  (22) 


|r-f0(s)|4  |  ds 


Similar  to  the  parallel-beam  case  described  above,  using  Eq.  (15),  one  can  conclude  that 


Cov{g(xc,  sa,  sb),  g(x'c,  sa,  Sh)}  «  c(xc)  6(xc  -  x' ), 


(23) 


where 


c(xc) 


A2{ud) 

-ro{s)  I 


4 


dr0(s) 

ds 


2 

<72(ud,s). 


(24) 


3.1.3  Noise  properties  of  the  differentiation/backprojection  for  cone-beam  data 

In  the  cone-beam  case,  let  P(itd,  Vd,  s)  denote  the  projection  at  view  s  on  a  detector  bin  specified  by  ( Ud ,  Vd). 
In  the  so-called  “well  oriented"  rotation-coordinate  system  [5],  unit  vector  eu(s)  is  parallel  to  and  unit  vectors 
e„  (s)  and  ew  ( s )  are  orthogonal,  to  the  tangential  direction  di^(s)  0f  the  source  trajectory.  Let  Ud  and  Vd  denote 
the  coordinates  along  e„(s)  and  ev(s).  It  can  be  shown  [5]  that  the  backprojection  image  depends  only  upon  the 
data  derivative  along  Ud .  Therefore,  the  reconstruction  formula  for  the  cone-beam  backprojection  image  can  be 
obtained  from  that  for  the  fan-beam  backprojection  image  in  Eq.  (18)  by  simply  replacing  P(iid,  s)  and  A(v,d ) 
with  P(ud,  i’d,  s)  and  A(ud,  Vd)~  respectively,  in  Eq.  (19).  Subsequently,  one  can  show  that  the  covariance  of  the 
cone-beam  backprojection  image  g(xc,  sa,  sb)  also  satisfies  Eqs.  (23)  and  (24). 


3.1.4  Estimation  of  parameters  a  and  to  in  discrete  form 

The  parameter  a  is  introduced  to  account  for  the  interpolation  effect  of  the  discrete  data-derivative  on  the  chord- 
image  variance.  We  consider  a  two-point  derivative,  which  was  used  in  our  numerical  studies.  Let  P,  denote  the 
discrete  data,  where  *  =  1,2, ...,  /,  and  /  indicates  the  total  number  of  data  points.  We  assume  that  data  P,  are 
uncorrelated  and  use  Var{P, }  to  denote  their  variances.  The  two-point  data  derivative  is  defined  as 

Pi  =  ^[Pi+i-Pi-i].  (25) 

Therefore,  the  variance  of  the  discrete  data  derivative  P,  can  be  written  as 

VMP.)  .  IVarm+.HVari,,-.}  (26) 

When  data  variances  are  identical,  Eq.  (26)  becomes 

Var{Pj}  -  ivar{Pj.  (27) 

Therefore,  in  our  studies,  we  select  a  =  which  is  the  coefficient  in  Eq.  (27). 

The  parameter  u>  was  introduced  to  account  for  the  interpolation  effect  of  discrete  backprojection  on  the  chord- 
image  variances.  The  estimated  value  of  u>  depends  obviously  upon  the  specific  interpolation  scheme  used  in 
the  discrete  backprojection.  We  illustrate  below  our  estimation  of  oj  when  a  two-point  interpolation  is  used  for 
the  discrete  backprojection  in  the  parallel-  and  fan-beam  cases.  One  can  readily  obtain  estimates  of  u>  when 
other  interpolation  schemes  are  used.  At  a  backprojection  view  s,  we  use  Pj  to  denote  the  discrete  weighed-data 
derivatives.  For  a  given  Ud  satisfying  i,  <  Ud  <  (i  +  1),  we  express  the  interpolated  weighted-data  derivative  as 

P«rf  =  (l  — 7)P'+7P'+i,  (28) 

where  7  =  Ud  —  *.  Furthermore,  we  can  write  the  variance  of  P'Ud  as 


Var{P;d}«[(l-7)2+72]Var{P'} 


(29) 


For  the  seek  of  simplifying  the  estimation  of  u>,  we  have  ignored  the  correlation  between  P-  and  P'+1  and  assumed 
that  Var{P'}  «  Var{Pj+1}.  We  select  to  as  the  average  over  all  of  the  possible  7’s,  which  can  be  computed  as 

u  =  Jq  K1  ~  T')2  +  72]^7  =  (30) 

Finally,  substitution  of  a  =  1  and  oj  =  |  into  Eqs.  (17)  and  (24),  we  obtain  the  variances  c(xc)  of  the  backpro- 
jection  images  on  the  chords  for  the  parallel-beam  and  fan-beam  projections,  respectively,  as 


c(xc) 


ds  a2(ud,s), 


c(xc) 


A2{ud) 

|r-f0(s)|4 


dr0(s ) 
ds 


2 

a2(ud,s). 


(31) 

(32) 


3.1.5  Noise  property  of  the  weighted  Hilbert  transform  over  a  finite  interval 


The  weighted  Hilbert  transform  constitutes  an  important  step  in  the  chord-based  BPF,  MDFBP,  and  FBP  algo¬ 
rithms.  Consequently,  the  noise  properties  of  these  algorithms  depend  upon  that  of  the  weighted  Hilbert  trans¬ 
form,  which  we  study  below.  Let  f (xc .  sa,  sb)  denote  the  weighted  Hilbert  transform  of  the  backprojection  image 

g(xc,sa,sb ): 


f  [Xci  Sai  Sb)  — 


dx' 


b{xc) 


Xc  ~  X' 


■b{x'c)g{x'c,sa,sb). 


(33) 


We  assume  that  g(x'c,  s0,  sb)  is  band-limited  to  .  Therefore,  the  Hilbert  transform  kernel  -j-  can  be  replaced  by 


h{xc)  =  -7 rj 


dv  sgn[z7]e27UI'Xc 


2sin2(7rt/mxc) 

xc 


(34) 


In  the  presence  of  noise,  the  weighted  Hilbert  transform  i(xc,  sa,  sb)  should  be  interpreted  as  a  stochastic  process, 
which  is  denoted  in  boldface.  The  variance  of  f(xc,  sa,  sb )  can  be  written  as 

1  CXB  CXb 

Var{f(atc,s0,sb)}  =  — — -  /  dx'c  /  dx"c  b(x'c)  b(x")  (35) 

o2(xc)  JXA  JXA 

x  h(xc  -  x’c)  h{xc  -  x")  Cov{g(a:^  sa,  sb),  g(x",  sa ,  sb)}. 


As  Eqs.  (16)  and  (23)  show,  the  backprojection  image  g(xc,  sa,  sb)  can  be  treated  as  an  approximated  uncorrelated 
stochastic  process.  Using  the  result  in  Eq.  (16)  or  (23),  we  can  write  Eq.  (35)  as 

-  1  fXB 

Var{f(xc,sa,sh)}  =  -rj- — -  dxc  c(x'c)  h2{xc  -  x'c)  b2(x'c).  (36) 

°2{xc)  JXA 

In  our  numerical  studies  in  Secs.  3.2.1  and  3.3.1,  we  have  used  where  Ac  denote  the  sample  interval 

Of  S(xc,  Sa,  Sb)- 


3.1.6  Noise  properties  of  chord  images 

Using  Eq.  (4),  one  can  write  the  variance  of  the  reconstructed  chord  image  as 


Var{fc(a;c,  sa,  sb)}  «  Var{f(a:c,  sa,  sb)} 


(37) 


+ 


Var{P(urf0,  vdQ,  sa)| 
47r262(a;c) 


\/((  —  xB){xA  +  0  y/JJ  +  xa){1  +  xb) 


~\  2 


l  +  Xc 


l  —  X( 


a  b  c 

Figure  1 :  (a)  Phantom  in  the  numerical  studies,  (b)  Parallel-beam  configuration,  (c)  Fan-beam  configuration.  The 
solid  line  segment  with  the  endpoints  xa  and  x b  represents  the  reconstruction  segment.  The  thick  segment  with 
endpoints  xci  and  xC2  indicate  the  support  segment.  It  should  be  noted  that  [xa,  xb]  2  [xci,xc2\-  The  rectangular 
ROI  is  decomposed  into  a  set  of  (dashed)  line  segments. 


Circular  parallel-beam 
scan 

Circular  fan-beam 

scan 

Helical  cone-beam 

scan 

Number  of  views  per  turn 

1024 

512 

300 

Angular  range 

r  71  711 

L  2  ’  2  J 

[-7T,  ir\ 

[— 7T,  7t] 

Detector  bins 

256 

256 

128  x  128 

Detector  bin  size  (mm) 

1.1 

1.6 

3. 1x3.1 

SDD  (mm) 

- 

200 

200 

SID  (mm) 

- 

200 

200 

Helical  pitch  (mm) 

- 

- 

90 

Table  1:  Parameters  for  circular  parallel-beam,  circular  fan-beam,  and  helical  cone -beam  scanning  configurations 
which  are  used  in  our  simulation  studies.  SDD  is  the  source-to-detector  distance  and  SID  denotes  the  source-to- 
isocenter  distance. 


Substituting  Eq.  (36)  into  Eq.  (37)  yields 

1  fXB 

Var{fc(xc,  sa,  Sb)}  «  757 — r/  dx'c  c(x'c)  h2(xc  -  x’c)  b2{x’c) 

b  (xc)  Jxa 

Var{P(tfrf0,  Vdo,  sa)}  \/{l-  xB){xA  +  l 
4:Tr2b2(xc)  l  —  xc 

which  provides  a  formula  for  computing  the  chord-image  variance. 

3.2  Numerical  studies  of  noise  properties  in  parallel-beam  reconstruction 

Using  the  parallel-beam  configuration  in  Fig.  lb  and  the  parameters  given  in  Table  1,  we  calculated  noiseless 
projections  for  the  numerical  phantom  in  Fig.  la.  We  have  used  an  object-independent  Gaussian  noise  model  and 
an  object-dependent  Poisson-noise  model  in  the  numerical  studies.  For  each  noise  model,  we  generated  10,000 
sets  of  noisy  data  by  using  noiseless  data  as  the  means.  The  standard  deviation  cto  of  Gaussian  noise  used  is  1 .6% 
of  the  maximum  value  in  the  noiseless  data,  whereas  the  standard  deviation  for  the  Poisson  noise  is  the  noiseless 
data  scaled  to  yield  a  total  photon  count  of  5  x  105  for  each  view.  We  investigated  four  reconstruction  segments 
with  different  lengths  Lab  =  \xb  —  xA |:  7.8  cm,  10.0  cm,  14.1  cm,  and  20.0  cm,  all  of  which  are  located  at 
x  =  4.06  cm.  It  can  be  observed  in  Fig.  lb  that  the  length  of  the  support  segment,  5.5  cm  in  length,  is  shorter  than 
the  four  reconstruction  segments  considered.  Therefore,  the  image  on  this  chord  can  be  reconstructed  exactly  by 
use  of  data  determined  by  these  reconstruction  segments.  One  can  also  conclude  from  Fig.  lb  that  the  minimum 
data  required  by  the  first  three  reconstruction  segments,  which  are  shorter  than  the  maximum  dimension  (about 
15.6  cm)  of  the  object  support,  are  truncated. 


(38) 


0  yj  (l  +  xa){1  +  xb) 


l  T  xc 


3.2.1  Noise  properties  in  reconstruction  from  truncated  parallel-beam  data 

From  the  10,000  sets  of  data  containing  Gaussian  noise,  we  used  Eqs.  (4),  (5),  and  (6)  to  reconstruct  10,000 
noisy  fc(xc,  sa,  St),  f  (xc,  sa,  s&),  and  g(xc,  sa,  S&),  respectively,  on  the  four  reconstruction  segments  described 
above.  Based  upon  these  noisy  reconstructions,  we  subsequently  computed  their  corresponding  empirical  vari¬ 
ances,  which  are  shown  in  the  upper  row  of  Fig.  2.  We  compare  the  empirical  results  to  the  analytical  results 
obtained  by  use  of  Eqs.  (16),  (36)  and  (38).  The  function  c(xc)  is  determined  by  using  cr(ud,  s)  =  ctq  in  Eq.  (31), 
where  op  is  1 .6%  of  the  maximum  value  in  noiseless  data.  The  analytical  results  are  displayed  in  the  lower  row 
of  Fig.  2.  Similarly,  using  Eqs.  (4),  (5),  and  (6),  we  reconstructed  10,000  sets  of  noisy  g(xc,  sa,  Sb),  f(xc,  sa,  Sb), 
and  ic(xc,  sa,  Sb )  on  the  four  segments  from  10,000  sets  of  data  containing  Poisson  noise.  The  computed  empir¬ 
ical  variances  from  these  noisy  images  are  displayed  in  the  upper  row  of  Fig.  3.  Using  the  noiseless  data  as  the 
Poisson-noise  variance  cr2(ud ,  s)  in  Eq.  (17),  one  can  readily  determine  c(xc);  and  using  the  determined  c(xc)  in 
Eqs.  (16),  (36)  and  (38),  one  can  compute  analytical  image  variances,  which  are  displayed  in  the  lower  row  of 
Fig.  3.  The  results  show  that  the  analytical  and  empirical  results  agree  well  with  each  other,  suggesting  that  Eq. 
(38)  provides  an  adequate  estimation  of  the  chord-image  variance. 

It  can  also  be  observed  in  Figs.  2c  and  3c  that,  the  shorter  the  reconstruction  segment,  the  higher  the  chord- 
image  variances.  This  is  only  because  the  second  term  in  (38)  increases  as  Lab  (i.e.,  (xa  —  xc)(xc  —  Xb)) 
decreases.  However,  the  difference  of  the  chord-image  variances  in  the  central  part  of  the  support  segment  is  quite 
small  among  these  reconstruction  segments.  The  implication  of  this  result  is  that  there  may  be  a  significant  gain 
in  terms  of  dose  reduction  by  using  a  short  reconstruction  segment,  because  data  required  to  reconstruct  an  image 
on  this  reconstruction  segment  is  less  than  that  required  by  using  a  longer  reconstruction  segment,  thus  resulting 
in  a  reduced  illumination  coverage  to  the  object.  For  similar  x-ray  intensities,  which  is  directly  related  to  the 
data-noise  level,  the  reconstruction  using  a  short  reconstruction  segment  appears  to  yield  image  variance  within 
the  support  segment  that  is  comparable  to  that  obtained  with  a  longer  reconstruction  segment. 

We  have  also  performed  numerical  studies  of  the  noise  properties  of  the  reconstructed  ROI-images  by  use  of  the 
BPF  and  MDFBP  algorithms  from  truncated  data.  Using  the  numerical  phantom  in  Fig.  1  and  each  of  Gaussian- 
and  Poisson-noise  models  described  above,  we  generated  500  noisy,  truncated  data  sets  for  image  reconstruction 
on  reconstruction  segments  of  a  length  Lab=^0.0  cm,  as  shown  in  Fig.  la,  which  completely  cover  the  ROI. 
We  subsequently  reconstructed  500  noisy  images  by  using  the  BPF  and  MDFBP  algorithms.  We  display  in  Fig. 
4  noisy  ROI  images  reconstructed  using  the  BPF  and  MDFBP  algorithms  from  data  containing  Gaussian  noise 
(upper  row)  and  Poisson  noise  (lower  row). 

Using  the  reconstructed  500  sets  of  Gaussian-noise  images  and  500  sets  of  Poisson-noise  images,  we  computed 
empirical  variances  of  the  ROI  images,  which  are  shown  in  the  upper  row  and  lower  row  of  Fig.  5,  respectively. 
We  display  in  the  third  column  of  Fig.  5  the  variance  profiles  on  the  dashed  lines  indicated  in  the  variance  images. 
Results  in  Fig.  5  support  the  conclusion  that  both  BPF  and  MDFBP  algorithms  yield  images  with  comparable 
variance  levels. 

3.2.2  Noise  properties  in  reconstruction  from  non-truncated  parallel-beam  data 

As  discussed  above,  the  FBP  algorithm  cannot  reconstruct  exactly  images  from  truncated  data.  Therefore,  we 
study  below  the  noise  properties  of  the  FBP  algorithm  from  parallel-beam  data  without  truncations.  For  the 
purpose  of  comparison,  we  have  also  included  reconstruction  results  of  the  BPF  and  MDFBP  algorithms  from 
the  same  non-truncated  data.  Using  the  numerical  phantom  in  Fig.  1  and  each  of  the  Gaussian-  and  Poisson- 
noise  models,  we  generated  500  noisy  data  sets  from  which  500  noisy  images  were  obtained  by  use  of  each 
of  the  BPF,  MDFBP,  and  FBP  algorithms.  Using  these  noisy  images,  we  computed  empirical  variance  images, 
which  are  shown  in  the  upper  and  lower  rows  of  Fig.  6,  respectively,  for  the  Gaussian-  and  Poisson-noise  models. 
We  also  display  in  Fig.  7  the  variance  profiles  on  the  dashed  lines  (i.e.,  on  a  chord)  indicated  on  the  variance 
images  in  Fig.  6.  The  profile  results  were  obtained  by  use  of  the  BPF  (solid),  MDFBP  (dashed),  and  FBP  (dotted) 
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Figure  2:  Empirical  (top  row)  and  analytical  (middle  row)  variances  of  (a)  g(xc,  sa,  Sb),  (b)  f(xc,sa,Sb),  and 
(c)  fc(xc,  sa,  Sb)  obtained  on  four  reconstruction  segments  from  parallel-beam  data  containing  Gaussian  noise. 
The  difference  between  empirical  variances  and  analytical  variances  is  also  shown  in  the  bottom  row,  which 
demonstrates  the  analytical  variances  agree  well  with  the  empirical  variances.  The  lengths  of  these  segments  are 
indicated  in  the  box  in  upper-right  corners  of  the  plots. 
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Figure  3:  Empirical  (top  row)  and  analytical  (middle  row)  variances  of  (a)  g (xc,  sa,  Sb),  (b)  f(xc,sa,Sb),  and 
(c)  f c(xc,  sa,  Sb)  obtained  on  four  reconstruction  segments  from  parallel-beam  data  containing  Poisson  noise. 
The  difference  between  empirical  variances  and  analytical  variances  is  also  shown  in  the  bottom  row,  which 
demonstrates  the  analytical  variances  agree  well  with  the  empirical  variances.  The  lengths  of  these  segments  are 
indicated  in  the  box  in  upper-right  corners  of  the  plots. 
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Figure  4:  ROI-images  obtained  by  use  of  the  BPF  algorithm  (a)  and  MDFBP  algorithm  (b)  from  truncated  parallel- 
beam  data  containing  Gaussian  noise  (upper  row)  and  Poisson  noise  (lower  row). 


Figure  5:  Empirical  ROI-image  variances  obtained  by  use  of  the  BPF  (a)  and  MDFBP  (b)  algorithms  from  trun¬ 
cated  parallel-beam  data  containing  Gaussian  noise  (upper  row)  and  Poisson  noise  (lower  row),  respectively,  (c) 
Variance  profiles  on  the  dashed  lines  indicated  in  columns  (a)  and  (b)  obtained  with  the  BPF  (solid)  and  MDFBP 
(dashed)  algorithms,  respectively.  For  the  purpose  of  displaying  the  details  in  the  central  (i.e.,  low  variance)  re¬ 
gions,  we  have  applied  a  logarithmic  scale  to  the  variance  images.  The  display  windows  are  [-0.86,  0.50]  and 
[-1.28,  0.20]  for  Gaussian  noise  and  Poisson  noise,  respectively. 


a 


b 


c 


Figure  6:  Empirical  variance  images  obtained  by  use  of  the  BPF  (a),  MDFBP  (b),  and  FBP  (c)  algorithms  from 
non-truncated  parallel-beam  data  containing  Gaussian  noise  (upper  row)  and  Poisson  noise  (lower  row).  For  the 
purpose  of  displaying  the  details  in  the  central  (i.e.,  low  variance)  regions,  we  have  applied  a  logarithmic  scale  to 
the  variance  images.  The  display  windows  are  [-1.3,  0.30]  and  [-1.3,  0.40]  for  Gaussian  noise  and  Poisson  noise, 
respectively. 

algorithms  for  the  Gaussian-noise  model  (a)  and  Poisson-noise  model  (b).  It  can  be  observed  that  image  variances 
obtained  with  the  three  algorithms  are  similar  and  that  the  only  difference  comes  at  the  extreme  ends  of  the  shown 
reconstruction  segments.  The  BPF  and  MDFBP  algorithms  show  a  significant  increase  in  the  image  variance  at 
both  ends  of  the  profile.  The  reason  for  this  is  that  the  reconstruction  segment  was  taken  to  be  the  width  of  the 
image  array,  and  the  pre-factor  for  the  finite  Hilbert  transform  in  Eqs.  (4)  and  (7)  has  a  singularity  at  the  ends  of  the 
reconstruction  segment.  In  practical  situations  this  pre-factor  is  of  little  consequence  because  the  reconstruction 
segment  can  be  selected  larger  to  avoid  the  singular  behavior;  furthermore,  because  the  singularity  goes  as  the  -1/2 
power,  its  effect  is  evident  only  very  close  to  the  endpoints  of  the  reconstruction  segment. 

3.3  Numerical  studies  of  noise  properties  in  fan-beam  reconstruction 

Using  the  fan-beam  configuration  in  Fig.  lc  and  the  parameters  listed  in  Table  1,  we  calculated  fan-beam,  noiseless 
data  for  the  numerical  phantom  in  Fig.  la.  We  have  also  used  an  object-independent  Gaussian-noise  model  and  an 
object-dependent  Poisson-noise  model  in  this  numerical  study.  The  standard  deviation  op  of  Gaussian  noise  used 
is  2.3%  of  the  maximum  value  in  noiseless  fan-beam  data,  whereas  the  standard  deviation  for  the  Poisson  noise  is 
the  noiseless  data  scaled  to  yield  a  total  photon  count  of  5  x  105  for  each  view.  For  each  noise  model,  10,000  sets  of 
noisy  data  were  generated  by  use  of  the  corresponding  noiseless  data  as  the  means.  We  investigated  reconstruction 
segments  of  four  different  lengths  Lab  =7.8  cm,  10.0  cm,  14.1  cm,  and  20.0  cm.  All  of  the  segments  are  located 
at  x  =  4.06  cm.  It  can  be  observed  in  Fig.  lc  that  the  length  of  the  support  segment  is  5.5  cm,  which  is  shorter 
than  the  four  reconstruction  segments.  Therefore,  the  image  on  this  chord  can  be  reconstructed  exactly  by  use 
of  data  determined  by  these  reconstruction  segments.  One  can  also  conclude  from  Fig.  lc  that  data  determined 
by  the  first  three  reconstruction  segments,  which  are  shorter  than  the  maximum  dimension  (about  15.6  cm)  of  the 
object  support,  are  truncated. 
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Figure  7:  Variance  profiles  on  the  dashed  lines  indicated  in  variances  images  shown  in  Fig.  6.  They  were  obtained 
for  the  Gaussian  (a)  and  Poisson  (b)  noise  models  by  use  of  the  BPF  (solid),  MDFBP  (dashed),  and  FBP  (dotted) 
algorithms,  respectively. 

3.3.1  Noise  properties  in  reconstruction  from  truncated  fan-beam  data 

From  the  10,000  sets  of  data  containing  Gaussian  noise,  we  used  Eqs.  (4),  (5),  and  (6)  to  reconstruct  10,000 
noisy  fc(xc,  sa,  St),  f(xc,  sa,  Sb ),  and  g(xc,  sa,  s&),  respectively,  on  the  four  reconstruction  segments  described 
above.  Based  upon  these  noisy  reconstructions,  we  subsequently  computed  their  corresponding  empirical  vari¬ 
ances,  which  are  shown  in  the  upper  row  of  Fig.  8.  As  for  the  analytic  variance,  one  can  determine  c(xc)  by  using 
<j(ud ,  s)  =  (7 o  in  Eq.  (32),  where  op  is  2.3%  of  the  maximum  value  in  noiseless  fan-beam  data.  Using  c(xc)  in 
Eqs.  (23),  (36)  and  (38),  we  computed  analytically  image  variances,  which  are  displayed  in  the  lower  row  of  Fig. 
8.  Similarly,  using  Eqs.  (4),  (5),  and  (6),  we  reconstructed  10,000  sets  of  noisy  fc(xc,  sa,  St),  f(a;c,  sa,  Sb),  and 
g(xc,  sa ,  Sb)  on  the  four  segments  from  10,000  sets  of  fan-beam  data  containing  Poisson  noise.  The  computed 
empirical  variances  from  these  noisy  images  are  displayed  in  the  upper  row  of  Fig.  9.  Furthermore,  using  the 
noiseless  fan-beam  data  as  the  Poisson-noise  variance  a 2(ud,  s)  in  Eq.  (32),  one  can  readily  determine  c{xc). 
Using  the  determined  c(xc)  in  Eqs.  (23),  (36)  and  (38),  we  computed  analytically  image  variances,  which  are 
displayed  in  the  lower  row  of  Fig.  9.  It  can  be  observed  that  the  analytic  and  empirical  results  agree  well  with 
each  other,  suggesting  that  Eq.  (38)  provides  an  adequate  analytic  estimation  of  the  chord-image  variance  for  the 
fan-beam  case  as  well.  It  is  interesting  to  note  in  Figs.  8a  and  9a  that  the  variances  of  g(xc,  sa,  Sb )  is  spatially 
varying  on  the  chord.  Based  upon  Eq.  (32),  one  can  readily  conclude  that  this  spatial  variation  is  caused  by  the 
spatially  variant  factor  ^i^|4  . 

Again,  from  these  results,  observations  similar  to  those  for  the  parallel-beam  case  can  be  made  for  the  fan-beam 
case.  For  example,  as  Figs.  8c  and  9c  show,  the  shorter  the  reconstruction  segment,  the  higher  the  chord-image 
variances.  This  is  only  because  the  second  term  in  (38)  increases  as  Lab  (be.,  (xa  —  xc)[xc  —  xb))  decreases. 
The  implication  of  these  results  is  that  there  may  be  a  significant  gain  in  terms  of  dose  reduction  by  using  a  short 
reconstruction  segment,  because  data  required  to  reconstruct  an  image  on  this  reconstruction  segment  is  less  than 
that  required  by  using  a  longer  reconstruction  segment,  which  can  result  in  reduction  of  illumination  coverage  of 
the  object. 

We  have  also  performed  numerical  studies  of  the  noise  properties  of  the  reconstructed  ROI-images  by  use 
of  the  BPF  and  MDFBP  algorithms  from  truncated  data.  Using  the  numerical  phantom  in  Fig.  1  and  each  of 
the  Gaussian-  and  Poisson-noise  models  described  above,  we  generated  500  noisy  truncated  data  sets  for  image 
reconstruction  on  reconstruction  segments  of  a  length  Lab= 10.0  cm,  as  shown  in  Fig.  la,  which  cover  the  ROI 
completely.  We  subsequently  reconstruct  500  noisy  images  by  using  each  of  the  BPF  and  MDFBP  algorithms. 


pixel  position  (cm)  pixel  position  (cm)  pixel  position  (cm) 

a  b  c 


Figure  8:  Empirical  (top  row)  and  analytical  (middle  row)  variances  of  (a)  g(xc,  sa,  Sb),  (b)  f(xc,  sa,  Sb ),  and  (c) 
fc(xc,  sQ,  Sb)  obtained  on  four  reconstruction  segments  from  fan-beam  data  containing  Gaussian  noise.  The  dif¬ 
ference  between  empirical  variances  and  analytical  variances  is  also  shown  in  the  bottom  row,  which  demonstrates 
the  analytical  variances  agree  well  with  the  empirical  variances.  The  lengths  of  these  segments  are  indicated  in 
the  box  in  upper-right  corners  of  the  plots. 
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Figure  9:  Empirical  (top  row)  and  analytical  (middle  row)  variances  of  (a)  g(xc,  sa,  Sb),  (b)  f(xc,  sa,  Sb ),  and  (c) 
fc(xc,  sQ,  Sb)  obtained  on  four  reconstruction  segments  from  fan-beam  data  containing  Poisson  noise.  The  differ¬ 
ence  between  empirical  variances  and  analytical  variances  is  also  shown  in  the  bottom  row,  which  demonstrates 
the  analytical  variances  agree  well  with  the  empirical  variances.  The  lengths  of  these  segments  are  indicated  in 
the  box  in  upper-right  corners  of  the  plots. 
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Figure  10:  Empirical  ROI-image  variances  obtained  by  use  of  the  BPF  (a)  and  MDFBP  (b)  algorithms  from 
truncated  fan-beam  data  containing  Gaussian  noise  (upper  row)  and  Poisson  noise  (lower  row),  respectively,  (c) 
Variance  profiles  on  the  dashed  lines  indicated  in  columns  (a)  and  (b)  obtained  with  the  BPF  (solid)  and  MDFBP 
(dashed)  algorithms,  respectively.  For  the  purpose  of  displaying  the  details  in  the  central  (i.e.,  low  variance) 
regions,  we  have  applied  a  logarithmic  scale  to  the  variance  images.  The  display  windows  are  [-0.76,0.29]  and 
[-0.75,-0.04]  for  Gaussian  noise  and  Poisson  noise,  respectively. 


Using  the  reconstructed  500  sets  of  Gaussian-noise  images  and  500  sets  of  Poisson-noise  images,  we  computed 
empirical  variance  images  within  the  ROI,  which  are  shown  in  the  upper  row  and  lower  row  of  Fig.  10,  respectively. 
Moreover,  we  display  in  the  third  column  of  Fig.  10  the  variance  profiles  on  the  dashed  lines  indicated  in  the 
variance  images.  Results  in  Fig.  10  support  the  conclusion  that  both  BPF  and  MDFBP  algorithms  yield  images 
with  comparable  variance  levels. 

3.3.2  Noise  properties  in  reconstruction  from  non-truncated  fan-beam  data 

The  FBP  algorithm  cannot  reconstruct  exactly  images  from  truncated  fan-beam  data.  Thus,  we  evaluate  below  the 
noise  properties  of  the  FBP  reconstruction  from  non-truncated  fan-beam  data.  For  the  purpose  of  comparison,  we 
have  also  included  reconstruction  results  of  the  BPF  and  MDFBP  algorithms  from  the  same  data  sets.  Using  the 
numerical  phantom  and  fan-beam  configuration  described  above,  we  generated  non-truncated  fan-beam  data  at 
512  views  uniformly  covering  2tt.  Using  the  noiseless  data  as  the  means,  we  generated  500  sets  of  data  containing 
Gaussian  noise  and  500  sets  of  data  containing  Poisson  noise.  The  standard  deviation  for  the  Gaussian  noise  is 
2.3%  of  the  maximum  value  of  the  noiseless  data,  whereas  the  standard  deviation  for  the  Poisson  noise  is  the 
noiseless  data  scaled  to  yield  a  total  photon  count  of  5  x  105  for  each  view.  For  a  given  chord  specified  by 
sa  and  Sb,  one  can  reconstruct  its  image  from  data  acquired  over  the  right-side  trajectory  (i.e.,  s  G  [sQ,  s&]),  as 
shown  in  Fig.  11a.  Conversely,  one  can  also  reconstruct  the  chord  image  from  data  acquired  with  both  right-side 
trajectory  (i.e.,  s  £  [sa,  s;,])  and  left-side  trajectory  (i.e.,  s  G  [sb,  sa]),  as  shown  in  Fig.  1  lb.  In  chord-based  image 
reconstruction,  we  decompose  image  area  into  chords  parallel  to  the  vertical  direction,  and  the  source  scans  from 
sa  to  Sb  and  then  from  st,  to  sa,  as  shown  in  Fig.  11. 

For  each  chord  in  the  set  covering  the  image  area,  we  hi  st  reconstructed  the  images  by  use  of  the  BPF,  MDFBP, 
and  FBP  algorithms  from  data  containing  Gaussian  noise  acquired  over  the  right-side  trajectory  specified  by 
s  G  [sa,Sb].  Subsequently,  we  computed  empirically  chord-image  variances  from  these  noisy  reconstructions. 
By  assembling  the  chord-image  variances,  we  obtain  the  variance  images,  which  are  shown  in  the  upper  row  of 


a 


b 


i 


Figure  11:  (a)  The  right-side  trajectory  (solid)  of  the  chord  (thick)  specified  by  sa  and  s&.  (b)  The  right-side  trajec¬ 
tory  (solid)  and  left-side  trajectory  (solid)  of  the  chord  (thick)  specified  by  sa  and  Sb-  The  scanning  configuration 
in  (b)  corresponds  to  a  full,  fan-beam  scan. 
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Figure  12:  Empirical  variance  images  obtained  by  use  of  the  BPF  (a),  MDFBP  (b),  and  FBP  (c)  algorithms  from 
non-truncated  fan-beam  data  containing  Gaussian  (upper  row)  and  Poisson  (lower  row)  noise.  For  the  purpose  of 
displaying  the  details  in  the  central  (i.e.,  low  variance)  regions,  we  have  applied  a  logarithmic  scale  to  the  variance 
images.  The  display  window  is  [-1.0,  0.65]  and  [-2.02,  0.29]  for  Gaussian  noise  and  Poisson  noise,  respectively. 


Fig.  12,  for  the  BPF,  MDFBP  and  FBP  algorithms,  respectively.  Similarly,  from  data  containing  Poisson  noise, 
we  obtained  the  image  variances,  which  are  displayed  in  lower  row  of  Fig.  12. 

We  show  in  column  one  of  Fig.  13  the  image  variances  on  a  chord  specified  by  sa  =  tt/2  and  Sb  = 
7t/2  obtained  from  data  containing  Gaussian  noise  (upper  row)  and  Poisson  noise  (lower  row),  respectively.  As 
already  seen  above,  the  variance  increases  as  the  position  along  the  chord  near  the  source  trajectory.  There  is 
little  difference  between  the  three  algorithms.  Furthermore,  these  variance  images  have  similar  properties:  the 
chords  on  the  right  part  have  higher  and  more  non-uniform  image  variance  than  those  on  the  left  part  in  the  image 
area.  In  column  2  of  Fig.  13,  we  show  the  profiles  on  the  middle  points  across  the  vertical  chords  (i.e.,  on  the 
middle  horizontal  lines  in  the  variance  images  shown  in  Fig.  12).  The  results  reveal  that  some  difference  of 
the  MDFBP  result  from  the  BPF  and  FBP  results  in  the  peripheral  region.  This  difference  may  be  attributed  to 
the  data  weighting  prior  to  the  backprojection  step,  which  differs  from  that  in  the  BPF  and  FBP  algorithms.  We 
should  point  out  that  this  difference  is  only  seen  in  the  extreme  periphery  of  the  imaging  area.  For  most  practical 
situations  these  three  algorithms  perform  virtually  identically  in  terms  of  image  variance. 

We  investigate  further  the  general  trend  of  the  variance  decreasing  for  chords  on  the  left  of  the  variance  image. 
The  first  impression  of  this  behavior  is  that  this  trend  is  obvious,  because  the  scanning  trajectory  is  on  the  right 


Figure  13:  Variance  profiles  along  the  central  vertical  (a)  and  horizontal  (b)  lines  in  variance  images  shown 
in  Fig.  12,  obtained  by  the  BPF  (solid),  MDFBP  (dashed),  and  FBP  (dotted)  algorithms  from  data  containing 
Gaussian  noise  (upper  row)  and  Poisson  noise  (lower  row). 

side,  chords  on  the  left  of  the  variance  image  are  reconstructed  with  a  longer  scanning  trajectory.  It  appears  that 
more  data  are  used  in  reconstructing  chords  covering  the  left  part  of  the  variance  image.  This  explanation  is, 
however,  incorrect.  First,  there  is  a  slight  upturn  in  the  variance  for  the  chords  on  the  extreme  left  of  the  variance 
image,  which  runs  counter  to  this  trend.  Second,  it  can  be  demonstrated  that  the  amount  of  data  going  into  the 
chord  reconstruction  does  not  necessarily  increase  as  the  scanning  trajectory  increases.  Based  upon  Eq.  (32),  one 
can  conclude  that  the  true  cause  of  the  variance  behavior  is  spatially  dependent  weighting  factor,  l/|r  —  ro(A)|,  in 
the  BPF,  MDFBP,  and  FBP  algorithms  [19]. 

For  a  given  chord  specified  by  sa  and  Sb,  when  full  scan  data  are  available,  one  can  reconstruct  two  chord 
images  by  use  of  data  acquired  with  the  right-side  and  left-side  trajectories,  as  shown  in  Fig.  1 1  and  then  obtain 
a  final  chord  image  by  averaging  the  two  chord  images.  We  show,  in  Fig.  14,  the  variance  images  of  the  full  scan 
with  accompanying  profiles  in  Fig.  15. 

3.4  Numerical  studies  of  noise  properties  in  cone-beam  reconstruction 

The  BPF,  MDFBP,  and  FBP  algorithms  can  yield  exact  image  reconstruction  on  a  chord  specified  by  sa  and 
Sb  as  long  as  the  support  segment  on  the  chord  is  illuminated  by  the  x-ray  beam  at  the  projection  views  s  & 
[aa,  Sb],  because  these  algorithms  require  data  only  on  the  fan-beam  projections  of  the  support  segment.  From  the 
perspective  of  the  chord-based  algorithms,  the  reconstruction  of  a  chord  image  from  cone-beam  data  is  similar  to 
that  of  a  chord  image  from  fan-beam  data.  In  the  fan-beam  case,  the  orientation  of  the  fan-beam  planes  at  different 
views  remain  unchanged,  whereas,  in  the  cone-beam  case  with  a  non-planar  trajectory,  the  orientation  of  the  fan- 
beam-illumination  planes  generally  varies  from  view  to  view.  As  discussed  in  Sec.  3.1.3,  the  noise  properties 
of  differentiation,  backprojection  and  filtration  in  the  cone-beam  case  are  similar  to  that  in  the  fan-beam  case. 
Therefore,  we  include  below  only  the  study  results  on  the  noise  properties  of  the  final  chord-images  reconstructed 
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Figure  14:  Variance  images  obtained  by  use  of  the  BPF  (a),  MDFBP  (b),  and  FBP  (c)  algorithms  from  full-scan 
fan-beam  data  containing  Gaussian  (upper  row)  and  Poisson  (lower  row)  noise.  For  the  purpose  of  displaying  the 
details  in  the  central  (i.e.,  low  variance)  regions,  we  have  applied  a  logarithmic  scale  to  the  variance  images.  The 
display  windows  are  [-1.19,  0.31]  and  [-1.18,  0.20]  for  Gaussian  noise  and  Poisson  noise,  respectively. 


Figure  15:  Variance  profiles  along  the  central  vertical  (a)  and  horizontal  (b)  lines  in  variance  images,  which 
are  shown  in  Fig.  14,  obtained  by  the  BPF  (solid),  MDFBP  (dashed),  and  FBP  (dotted)  algorithms  from  data 
containing  Gaussian  noise  (upper  row)  and  Poisson  noise  (lower  row). 
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Figure  16:  (a)  7r-line,  37r-line,  and  57r-line  segments  in  a  helical  scan,  (b)  Surfaces  generated  in  the  imaging 
volume  by  concatenating  7r-line  segments  specified  by  sa  =  —  7r  and  Sb  £  [— 0.57T,  0.57t]. 

from  cone-beam  data. 

3.4.1  Helical  cone-beam  configuration 

In  our  investigation  of  the  noise  properties  of  image  reconstruction  from  cone-beam  data,  we  consider  the  helical 
trajectory,  which  is  the  most  widely  used  in  clinical  and  industrial  CT.  For  a  helical  scan,  the  source  trajectory  is 
described  mathematically  as  f*o(.s)  =  ( R  cos  s,  R  sin  s,  j^s),  where  R  is  the  source  to  center-of-rotation  distance, 
and  h  indicates  the  helical  pitch  length.  For  a  chord  specified  by  sa  and  S&,  if  (n  —  l)7r  <  |s&  —  sa|  <  (n  +  l)7r, 
where  n  is  a  positive  odd  integer,  the  chord  is  also  referred  to  as  an  n7r-line  segment  [20,  21],  as  shown  in  Fig. 
16a.  In  particular,  when  n  =  1  and  thus  0  <  |s&  —  sa|  <  2-7T,  the  chord  is  referred  to  as  a  7r-line  segment  [2,  11]. 
In  this  work  we  consider  image  reconstruction  only  on  7r-line  segments  for  the  reason  that  the  imaging  volume 
enclosed  by  the  helix  can  be  filled  uniquely  and  completely  by  7r-line  segments  [11,  12],  Thus,  7r-line  segments 
can  be  used  to  form  3D  images  in  a  helical  cone-beam  scan. 

We  computed  the  noiseless  data  from  a  Shepp-Logan  phantom  by  use  of  the  configuration  parameters  in  Table 
1.  Using  the  noiseless  data  as  the  means,  we  subsequently  generated  500  sets  of  data  containing  Gaussian  noise 
and  500  sets  of  data  containing  Poisson  noise,  respectively.  The  standard  deviation  of  Gaussian  noise  is  chosen  to 
be  0.7%  of  the  maximum  value  of  the  noiseless  data,  whereas  the  standard  deviation  for  the  Poisson  noise  is  the 
noiseless  data  scaled  to  yield  a  total  count  of  5  x  105  for  each  view. 

3.4.2  Noise  properties  in  reconstruction  from  helical  cone-beam  data 

A  curved  surface  in  the  helix  volume  can  be  formed  by  a  set  of  7r-line  segments  for  which  we  fix  one  end-point 
at  sa  and  sweep  the  other  endpoint  over  a  range  Sb  £  [.smln.  smax].  We  show  in  Fig.  16a  a  curved  surface  obtained 
by  concatenating  a  set  of  7r-line  segments  specified  by  sa  =  —  7r  and  Sb  £  [— 0.57T,  0.57t].  Using  generated  noisy 
helical  cone-beam  data,  we  reconstructed  noisy  images  on  the  7r-line  surface  by  use  of  the  BPF,  MDFBP,  and  FBP 
algorithms.  From  these  noisy  images  we  subsequently  computed  empirical  image  variances  on  the  7r-line  surface. 
In  Fig.  17,  we  display  the  image  variances  obtained  with  BPF,  MDFBP  and  FBP  algorithms  from  data  containing 
Gaussian  noise  and  Poisson  noise. 

We  also  display  in  Figs.  18a  and  b  image  variances  on  the  7r-line  segment  specified  by  sa  =  —7 r  and  Sb  =  0  in 
the  surface,  obtained  from  data  containing  Gaussian  noise  and  Poisson  noise,  respectively.  The  image  variances 
show  similar  characteristics  to  that  of  fan-beam  image  variances  observed  in  Fig.  13.  Namely,  the  variance  image 
on  the  7r-line  surface  in  Fig.  17  has  a  structure  that  is  similar  to  the  right-side  scan  fan-beam  results  presented 
in  Sec.  3.3.1;  the  images  on  7r-line  segments  reconstructed  from  smaller  helix  segments  tend  to  have  higher  and 
more  non-uniform  variances.  The  similarity  with  the  fan-beam  case  is  not  surprising  because  the  geometrical 
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Figure  17:  Empirical  variance  images  on  the  7r-line  surface  shown  in  Fig.  16b  obtained  by  use  of  the  BPF  (a), 
MDFBP  (b),  and  FBP  (c)  algorithms  from  data  containing  Gaussian  noise  (upper  row)  and  Poisson  noise  (lower 
row),  respectively.  For  the  purpose  of  displaying  the  details  in  the  central  (i.e.,  low  variance)  regions,  we  have 
applied  a  logarithmic  scale  to  the  variance  images.  The  display  windows  are  [-1.80,  0.30]  and  [-1.80,  0.60]  for 
Gaussian  noise  and  Poisson  noise,  respectively. 

arrangement  of  the  7r-line  with  respect  to  its  scanning  trajectory  is  very  similar  to  the  relationship  between  the 
chords  and  corresponding  fan-beam  scanning  trajectory.  The  only  difference  is  that  there  is  an  out-of-plane  bend 
to  the  helix  segment. 

Regarding  the  non-uniform  shape  of  the  variance,  one  can  attribute  the  high  variance  in  the  image  periphery 
to  the  weighting  factors  multiplying  the  data  derivatives  before  backprojection.  As  the  algorithms  are  essentially 
the  same  for  chord-image  reconstruction  in  fan-  and  cone -beam  cases,  this  conclusion  should  come  as  no  surprise. 
In  the  2D  fan-beam  case,  the  variance  non-uniformity  and  level  was  reduced  by  equally  weighting  reconstructions 
for  both  left  and  right  side  scans  for  each  chord  of  the  scanning  circle.  For  the  helical  configuration,  it  is  clear  that 
in  a  typical  scan  there  will  be  some  overscan  for  nearly  all  the  chords  comprising  the  volume.  But  the  overscan 
part  of  the  trajectory  does  not  form  a  closed  loop  so  using  the  overscan  data  to  reduce  image  variance  is  not  as 
obvious  as  the  case  of  the  circular  scan.  Future  work  will  focus  on  how  to  utilize  the  overscan  data  for  non-closed 
trajectories  for  the  purpose  of  reducing  the  impact  of  data  noise  on  chord-base  ROI-image  reconstruction. 

4  Discussion 

In  this  article,  we  have  performed  analytic  and  numerical  investigation  of  the  noise  properties  of  chord-based 
image  reconstructions  from  parallel-,  fan-,  and  cone-beam  data.  One  of  the  main  points  of  the  investigation 
was  to  test  whether  or  not  the  reduced  illumination  in  designing  a  minimal  data  set  for  a  particular  ROI  leads  to  a 
significant  reduction  in  exposure.  The  idea  was  to  compare  the  statistical  properties  of  the  ROI  image  reconstructed 
from  noise  realizations  of  the  minimal  data  set  with  noise  realizations  of  the  full  data  set.  Similar  noise  levels 
were  used  in  both  data  sets,  which  are  equivalent  to  modeling  similar  incident  x-ray  beam  intensities.  Our  study 
indicates  that  the  resulting  image  variance  was  almost  the  same  for  images  reconstructed  from  both  data  sets. 
Thus,  the  minimal  data  set  for  ROI  reconstruction  leads  to  a  significant  overall  dose  reduction,  because  the  body 
is  exposed  to  lower  amount  of  ionizing  radiation  in  the  reduced  scan.  For  fan-beam  and  cone -beam  imaging, 
we  explored  the  noise  properties  of  the  extreme  periphery  of  the  imaging  region  by  investigating  large  fan-  and 
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Figure  18:  Variances  on  the  central  vertical  line  segment,  specified  by  sa  =  —7 r  and  Sb  =  0,  obtained  with  the 
BPF  (solid),  MDFBP  (dashed),  and  FBP  (dotted)  algorithms  from  data  containing  Gaussian  noise  (a)  and  Poisson 
noise  (b). 

cone-angles.  Image  variance  non-uniformity  was  found  to  be  caused  by  spatially  dependent  weighting  factors  in 
the  chord-based  reconstruction  algorithms.  This  work  represents  a  study  of  the  noise  properties  of  chord-based 
reconstruction  and  of  the  impact  of  physical  factors  on  ROI  imaging  in  fan-beam  and  cone-beam  CT.  In  seeking 
ways  to  reduce  the  impact  of  noise  in  volume  imaging,  we  will  investigate  schemes  to  incorporate  overscan 
data.  The  analysis  presented  in  this  work  can  directly  be  applied  to  chord-based  image  reconstruction  for  general 
trajectories.  Finally,  it  is  important  to  investigate  the  behavior  of  the  ROI-reconstruction  algorithms  when  other 
important  factors  are  included  in  the  data  model  such  as  x-ray  polychromaticity  and  non-linear  partial  volume 
averaging. 
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